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ABSTRACT 

If a large Quantum Computer (QC) existed today, what type of physical problems could we efficiently simulate on 
it that we could not simulate on a classical Turing machine? In this paper we argue that a QC could solve some 
relevant physical "questions" more efficiently. The existence of one-to-one mappings between different algebras 
of observables or between different Hilbert spaces allow us to represent and imitate any physical system by any 
other one (e.g., a bosonic system by a spin-1/2 system). We explain how these mappings can be performed 
showing quantum networks useful for the efficient evaluation of some physical properties, such as correlation 
functions and energy spectra. 
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1. INTRODUCTION 

Quantum simulation of physical systems on a QC has acquired importance during the last years since it is 
believed that QCs can simulate quantum physics problems more efficiently than their classical analogues 1 : The 
number of operations needed for deterministically solving a quantum many-body problem on a classical computer 
(CC) increases exponentially with the number of degrees of freedom of the system. 

In quantum mechanics, each physical system has associated a language of operators and an algebra realizing 
this language, and can be considered as a possible model of quantum computation. 2 As we discussed in a 
previous paper, 3 the existence of one-to-one mappings between different languages (e.g., the Jordan- Wigner 
transformation that maps fermionic operators onto spin-I/2 operators) and between quantum states of different 
Hilbert spaces, allows the quantum simulation of one physical system by any other one. For example, a liquid 
nuclear magnetic resonance QC (NMR) can simulate a system of 4 He atoms (hard-core bosons) because an 
isomorphic mapping between both algebras of observables exists. 

The existence of mappings between operators allows us to construct quantum network models from sets 
of elementary gates, to which we map the operators of our physical system. An important remark is that 
these mappings can be performed efficiently: we need a number of steps that scales polynomially with the 
system size. However, this fact alone is not sufficient to establish that any quantum problem can be solved 
efficiently. One needs to show that all steps involved in the simulation (i.e., preparation of the initial state, 
evolution, measurement, and measurement control) can be performed with polynomial complexity. For example, 
the number of different eigenvalues in the two-dimensional Hubbard model scales exponentially with the system 
size, so QC algorithms for obtaining its energy spectrum will also require a number of operations that scales 
exponentially with the system size. 3 

Typically, the degrees of freedom of the physical system over which we have quantum control constitute 
the model of computation. In this paper, we consider the simulation of any physical system by the standard 
model of quantum computation (spin-1/2 system), since this might be the language needed for the practical 
implementation of the quantum algorithms (e.g., NMR). Therefore, the complexity of the quantum algorithms is 
analyzed from the point of view of the number of resources (elementary gates) needed for their implementation 
in the language of the standard model. Had another model of computation being used, one should follow the 
same qualitative steps although the mappings and network structure would be different. 

The main purpose of this work is to show how to simulate any physical process and system using the least 
possible number of resources. We organized the paper in the following way: In section 2 we describe the 
standard model of quantum computation (spin-1/2 system). Section 3 shows the mappings between physical 
systems governed by a generalized Pauli's exclusion principle (fermions, etc.) and the standard model, giving 



examples of algorithms for the first two steps (preparation of the initial state and evolution) of the quantum 
simulation. In section 4 we develop similar steps for the simulation of quantum systems whose language has an 
infinite-dimensional representation, thus, there is no exclusion principle (e.g., canonical bosons). In section 5 we 
explain the measurement process used to extract information of some relevant and generic physical properties, 
such as correlation functions and energy spectra. We conclude with a discussion about efficiency and quantum 
errors (section 6), and a summary about the general statements (section 7). 

2. STANDARD MODEL OF QUANTUM COMPUTATION 

In the standard model of quantum computation, the fundamental unit is the qubit, represented by a two level 
quantum system |a) = a|0) + For a spin-1/2 particle, for example, the two "levels" are the two different 
orientations of the spin, |f) = |0) and |J,) = |1). In this model, the algebra assigned to a system of iV-qubits 
is built upon the Pauli spin-1/2 operators crj, (J 3 V and a{ acting on the j-th qubit (individual qubit). The 

N 

commutation relations for these operators satisfy an su(2)j algebra defined by ([i, v, A = x, y, z) 

i=i 

WL,vt] = 2iS jk£^\<?{, (1) 

where e^x is the totally anti-symmetric Levi-Civita symbol. Sometimes it is useful to write the commutation 
relations in terms of the raising and lowering spin-1/2 operators 

a± = ^— 

Any operation on a QC is represented by a unitary operator U that evolves some initial state (boot-up 
state) in a way that satisfies the time-dependent Schrodinger equation for some Hamiltonian H. Any unitary 
operation (evolution) U applied to a system of N qubits can be decomposed into either single qubit rotations 
Rn{ti) = e~ l ^ a ^ by an angle d about the /j, axis or two qubits Ising interactions R z j >z k — e wa ^ . This is 
an important result of quantum information, since with these operations one can perform universal quantum 
computation. It is important to mention that we could also perform universal quantum computation with single 
qubit rotations and C-NOT gates 4 or even with different control Hamiltonians. The crucial point is that we need 
to have quantum control over those elementary operations in the real physical system. 

In the following, we will write down our algorithms in terms of single qubit rotations and two qubits Ising 
interactions, since this is the language needed for the implementation of the algorithms, for example, in a liquid 
NMR QC. Again, had we used a different set of elementary gates our main results still hold but with modified 
quantum networks. 

As an example of such decompositions, we consider the unitary operator U(t) = e lHt , where H = aa],ala^, 
represents a time-independent Hamiltonian. After some simple calculations 2,3 we decompose U into elementary 
gates (one qubit rotations and two qubits interactions) in the following way 

U(t) = e ra<T ^ CT "* = e-^^e^^^e^^e^^^e-^^e-^^^e^^. (3) 

This decomposition is shown in Fig. 1, where the quantum network representation is displayed. In the same 
way, we could also decompose an operator U'(t) = e~ M ' T » cr *' T » t using similar steps, by replacing a l x <-> <r* in the 
right hand side of Eq. 3. 

3. SIMULATION OF FERMIONIC SYSTEMS 

As discussed in the Introduction, quantum simulations require simulations of systems with diverse degrees of 
freedom and particle statistics. Fermionic systems are governed by Pauli's exclusion principle, which implies that 
no more than one fermion can occupy the same quantum state at the same time. In this way, the Hilbert space 
of quantum states that represent a system of fermions in a solid is finite-dimensional (2 N for spinless fermions, 
where N is the number of sites or modes in the solid), and one could think in the existence of one-to-one mappings 



(2) 




Figure 1. Decomposition of the unitary operator U(t) = e taa :, ° into elementary single qubit rotations and two qubits 
interactions. Time t increases from left to right. 

between the fermionic and Pauli's spin-1/2 algebras. Similarly, any language which involves operators with a 
finite-dimensional representation (e.g., hard-core bosons, higher irreps of su(2), etc.) can be mapped onto the 
standard model language. 5 

In the second quantization representation, the (spinless) fermionic operators cj (c i ) are defined as the creation 
(annihilation) operators of a fermion in the i-th mode (i — N). Due to the Pauli's exclusion principle and 

the antisymmetric nature of the fermionic wave function under the permutation of two fermions, the fermionic 
algebra is given by the following commutation relations 

{c t ,Cj} = 0, {c\,Cj} = % (4) 

where {, } denotes the anticommutator. 

The Jordan- Wigner transformation 6 is the isomorphic mapping that allows the description of a fermionic 
system by the standard model 



c i - (n -°t J 



(5) 




(6) 

where a 1 ^ are the Pauli operators defined in section 2. One can easily verify that if the operators cr^ satisfy the 
su(2) commutation relations (Eq. 1), the operators c\ and c i obey Eqs. 4. 

We now need to show how to simulate a fermionic system by a QC. Just as for a simulation on a CC, the 
quantum simulation has three basic steps: the preparation of an initial state, the evolution of this state, and 
the measurement of a relevant physical property of the evolved state. We will now explain the first two steps, 
postponing the third until section 5. 

3.1. Preparation of the initial state 

In the most general case, any quantum state \ip) of N e fermions can be written as a linear combination of Slater 
determinants \<j) a ) 

L 



where 



a=l 

N e 



\4> a ) = n 4 i vac > ( 8 ) 

with the vacuum state |vac) defined as the state with no fermions. In the spin language, |vac) = |J.| • • • J.). 

We can easily prepare the states \4> a ) by noticing that the quantum gate, represented by the unitary operator 

U m = e^ c ^ m ) (9) 

when acting on the vacuum state, produces c^O) up to a phase factor. Making use of the Jordan- Wigner 
transformation (Eqs. 5, 6), we can write the operators U m in the spin language 

m-l 

»2<c n 

U m = e ^ . (10) 

The successive application of N e similar unitary operators will generate the state \(f) a ) up to an irrelevant global 
phase. 

L 

A detailed preparation of the fermionic state = ^ g a \4> a ) can be found in a previous work. 2 The basic 

idea is to use L extra (ancilla) qubits, then perform unitary evolutions controlled in the state of the ancillas, 
and finally perform a measurement of the z-component of the spin of the ancillas. In this way, the probability 
of successful preparation of \ip) is 1/L. (We need of the order of L trials before a successful preparation.) 

Another important case is the preparation of a Slater determinant in a different basis than the one given 
before 

N e 

|<^> = ;Q4 |vac>, (11) 

where the fermionic operators d\ 's are related to the operators c] through the following canonical transformation 

l f = e lM ^ (12) 

with ~~d[ = (d\, d\, ■ ■ ■ , djy), = (c{, c\, ■ ■ ■ , cjy), and M is an N x N Hcrmitian matrix. Making use of 
Thouless's theorem, 7 we observe that one Slater determinant evolves into the other, \(f>p) = U\(f> a ), where the 

unitary operator U = e~ l c M c can be written in spin operators using the Jordan- Wigner transformation 
and can be decomposed into elementary gates, 3 as described in section 2. Since the number of gates scales 
polynomially with the system size, the state \<j>p) can be efficiently prepared from the state \4> a ). 



3.2. Evolution of the initial state 



The second step in the quantum simulation is the evolution of the initial state. The unitary evolution operator 
of a time-independent Hamiltonian H is U(t) = e lHt . In general, H = K + V with K representing the kinetic 
energy and V the potential energy. Since we usually have [K, V] ^ 0, the decomposition of U(t), written in the 
spin language through the Jordan- Wigncr transformation (Eqs. 5,6), in terms of elementary gates (one qubit 
rotations and two qubits interactions), becomes complicated. To avoid this problem, we instead use a Trotter 
decomposition, so the evolution during a short period of time (At = t/Af with At — ► 0) is approximated. To 
order O(At) (first order Trotter breakup) 

TV 

U(t) = l[U(At), (13) 

9=1 

[/(At) = e iHAt = e i(K+V)At ^ e iKA te iVAt_ (w) 

The potential energy V is usually a sum of commuting diagonal terms, and the decomposition of e lVAt into 
elementary gates is straightforward. However, the kinetic energy K is usually a sum of noncommuting terms of 
the form c\cj + S-c i (bilinear fermionic operators), so we need again to perform a Trotter approximation of the 

operator e lKAt . As an example of such a decomposition, we consider a typical term e 1 ^ c i c 3 +c i c ^ At (i < j) ; when 
mapped onto the spin language gives 

e k = i+1 = e k =' +1 e k = i+1 . (15) 

The decomposition of each term on the right hand side of Eq. 15 into elementary gates was already described 
in previous work. 3 In section 2 and Fig. 1, we also showed an example of such a decomposition for i = 1 and 
j = 3. It is important to mention that the required number of elementary gates scales polynomially with the 
length \ j — i\. Notice that this step is not necessary for bosonic systems since no string of a\ operators is involved 
(see section 4). 

The accuracy of this method increases as At decreases, so we might require a large number of gates to perform 
the evolution with small errors. To overcome this problem, one could use Trotter approximations of higher order 
in At. 8 

3.3. Generalization: simulation of anyonic systems 

The concepts described in sections 3.1 and 3.2 can be easily generalized to other more general particle statistics, 
namely hard-core anyons. By "hard-core", we mean that only zero or one particle can occupy a single mode 
(Pauli's exclusion principle). 

The commutation relations between the anyonic creation and annihilation operators a\ and dj, are given by 

[a^a^e = [a\,a]]e = 0, 
Ka]]_ e = 5 ij (l-(e- ie + l)n j ) , (16) 
[nj,at] = Sija] , 

(i < j) where nj = djdj, [A, B]$ = AB — e l8 BA, with < 6 < 2ir defining the statistical angle. In particular, 
6 = 7r mod(27r) corresponds to canonical spinless fermions, while 8 = mod(27r) represents hard-core bosons. 

In order to simulate this problem with a QC made of qubits, we need to apply the following isomorphic and 
efficient mapping between algebras 



n 3 

where the Pauli operators where defined in section 2, and since they satisfy Eq. 1, the corresponding 
commutation relations for the anyonic operators (Eqs. 16) are satisfied, too. 

We can now proceed in the same way as in the fermionic case, writing our anyonic evolution operator in terms 
of single qubit rotations and two qubits interactions in the spin-1/2 language. As we already mentioned, anyon 
statistics have fermion and hard-core boson statistics as limiting cases. We now relax the hard-core condition on 
the bosons. 



„i0 



= nr 

l<j 

= 1(1 + 4 



A9 



-erJ a_ 



(17) 



4. SIMULATION OF BOSONIC SYSTEMS 



Quantum computation is based on the manipulation of quantum systems that possess finite number of degrees of 
freedom (e.g., qubits). From this point of view, the simulation of bosonic systems appears to be impossible, since 
the non existence of an exclusion principle implies that the Hilbert space used to represent bosonic quantum sates 
is infinite-dimensional; that is, there is no limit to the number of bosons that can occupy a given mode. However, 
sometimes we might be interested in simulating and studying properties such that the use of the whole Hilbert 
space is unnecessary, and only a finite sub-basis of states is sufficient. This is the case for physical systems with 
interactions given by the Hamiltonian 



N 



II a ' l 3 b i b 3 + fa 



(18) 



where the operators b\ (6J create (destroy) a boson at site i, and n, = b\b i is the number operator. The space 
dimension of the lattice is encoded in the parameters ctij and . Obviously, the total number of bosons Np in 
the system is conserved, and we restrict ourselves to work with a finite sub-basis of states, where the dimension 
depends on the value of Np. 

The respective bosonic commutation relations (in an infinite-dimensional Hilbert space) are 



[b i ,b j ]=0, [b i ,b]]=5 i: 



(19) 



However, in a finite basis of states represented by {\n\, n,2, • • • , tin) with m — 0, - ■ ■ ,Np}, where Np is the 
maximum number of bosons per site, the operators b\ can have the following matrix representation 



b\ = t 



<g> • • • <g> 1 < 



i th factor 



(20) 



where <g> indicates the usual tensorial product between matrices, and the (Np + 1) x (Np + 1) dimensional 
matrices 11 and w are 



(21) 



It is important to note that in this finite basis, the commutation relations of the bosons b\ differ from the 
standard bosonic ones (Eq. 19) 5 
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(22) 



and clearly {b]) Np+1 = 0. 

As we mentioned in the Introduction, our idea is to simulate any physical system in a QC made of qubits. 
For this purpose, we need to map the bosonic algebra into the spin-1/2 language. However, since Eqs. 22 imply 
that the linear span of the operators b\ and b i is not closed under the bracket (commutator), a direct mapping 
between the bosonic algebra and the spin-1/2 algebra (such as the case of the Jordan- Wigner transformation 
between the fermionic and spin-1/2 algebra) is not possible. Therefore, we could think in a one-to-one mapping 
between the bosonic and spin-1/2 quantum states, instead of an isomorphic mapping between algebras. Let us 
show a possible mapping of quantum states. 

We start by considering only the i-th site in the chain. Since this site can be occupied with at most Np 
bosons, it is possible to associate an Np + 1 qubits quantum state to each particle number state, in the following 

way 



|0)i <-> 


tollj.2 ' 


• • In p }% 






lotlj.2 ' 


■ ■ In p )% 




|2)i ~ 


I0I1T2 ' 


■ ■ lN P )i 


(23) 


Np), <- 


I0I1I2 ' 


■ ■ LVp)» 





where \n)i denotes a quantum state with n bosons in site i. Therefore, we need N(N P + 1) qubits for the 
simulation (where N is the number of sites) . In Fig. 2 we show an example of this mapping for a quantum state 
with 7 bosons in a chain of 5 sites. 

By definition (see Eqs. 20, 21) b\ \n)i — y/rT+l \n+ l)i, so the operator 

N P -1 

b\ = V^+lcj^a 7 ? 1 ' 1 , (24) 

n=0 

where the pair (n, i) indicates the qubit n that represents the i-th site, acts in the Np + 1 qubits states of Eqs. 
23 as b\\i ■ ■ ■ J. n _iTnJ.n+i ' ' ' In p )z = \/ n + 1 Ho • • • lnT™+iln+2 • • • J.jVj»)»- Then, its matrix representation in 
this basis is the same matrix representation of b\ in the basis of bosonic states. Similarly, the number operator 
can be written 

N P n.i , 1 

^ = E n ^^' ( 25 ) 

n=0 

and act as n^lo ■ ■ ■ J. n -iTnWi ' ' ' lN P )i = n || ■■ ■ lntn+iln+2 ■ ■ ■ In p )i- Notice that [b\ , J2n=o a z' z ] = °> 
which means that these operators conserve the total z-component of the spin and, thus, always keep states 
within the same subspace. 

We can now write down the Hamiltonian in Eq. 18 in the spin-1/2 algebra as 

JV 

H = ctij b\bj + faj Thifij, (26) 

where the operators b\ (bj) are given by Eq. 24 and n, bt Eq. 25. In this way, we are able to obtain physical 
properties of the bosonic system (such as the mean value of an observable, the mean value of the evolution 
operator, etc.) in a QC made of qubits. It is important to note that the type of Hamiltonian given by Eq. 18 is 
not the only one that can be simulatable using the described method. The only constraint is a fixed maximum 
number of bosons per site (or mode) . 
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Figure 2. Mapping of the bosonic state |0 Q ), of a chain of 5 sites and 7 bosons, into a spin-1/2 state (Eq. 23). 
4.1. Preparation of the initial state 

As in the fermionic case, the most general bosonic state of an N sites quantum system with Np bosons can be 
written as a linear combination of product states like 



\<t> a ) = mir(bir---(tf N r» iva C >, 



(27) 



N 



where K is a normalization factor, n, is the number of bosons at site i (J^ n, = iVp), and |vac) is the boson 

vacuum state (no particle state). Using the mapping described in Eq. 23, we can write the vacuum state in 
the spin language as |vac) = | Toil ■ ■ ■ In p )i ® • • • ® | Toil ■ ■ ■ In p )n and \(f> a ) = |J. ■ • • ] ni ■ ■ ■ In p )i ® • •• ® 
1 4-0 • • ■ Trijv ' ' ' In p ) n (see Fig. 2 for an example). Therefore, the preparation of \<p a ) in a QC made of qubits is 
an easy process: only N spins are flipped from the fully polarized state, where all spins are pointing down. 

L 

The preparation of a bosonic initial state of the form = ^ g a \<f> a ) is realized as in the fermionic case. 

a=l 

Again, we need to add L ancillas (extra qubits), perform controlled evolutions on their states, and finally perform 
a measurement of an spin component. 2 

4.2. Evolution of the initial state 

The basic idea is to use the first order Trotter approximation (see the fermionic case) to separate those terms 
of the Hamiltonian that belong to the kinetic energy K, from the ones that belong to the potential energy V 
(H = K + V, [K,V]^0),le., 

e lHAt ~ e lKAt e lVAt . (28) 




Figure 3. Decomposition of the unitary operator U(t) = es"""//" 1 into single qubit rotations and two qubits interac- 
tions. Time t increases from left to right. 

In general, K is a sum of non commuting terms of the form b.b, + bjb k , and we need to perform another first 
order Trotter approximation to decompose it into elementary gates (in the spin language). Then, a typical term 
e l ( b i b i+ b j b i) t wnen mapped onto the spin language (Eq. 24) gives 

ft ^T 1 , , , , , 

exp[- £ ^(n + \)(n< + \) [«<X l+M + ^+ 1 ' i )(a^ ^ +lj + >V£ +lj ) 

n,n'= 

+((7»'V™+ 1 » < - a^a^Xa^ia^ 1 * - <'^<' +1 ^")]], (29) 

where Np is the number of bosons. The terms in the exponent of Eq. 29 commute with each other, so the 
decomposition into elementary gates becomes straightforward. As an example (see Fig. 3), we consider a system 
of two sites with one boson. We need then 2(1 + 1) = 4 qubits for the simulation, and Eq. 24 implies that 
b\ = a ! 1 ^ 1 and b\ = a°l\^ 2 . Then, jWt+tyi)* becomes 

exp^'^VS'M' 2 ) x exp(|«^V^ 2 ^ 2 ) x exp(|^V^^ 2 ^ 2 ) x exp^V^^ 2 ) (30) 

v ™,,/ l '_0, 1_0. 2^1. 2\ „ 0,1_1,1_0,2_1,2>i w ^ _0,1„.1.1„.0.2„.1,2\ w JOA ^1-1 „0,2„1, 2\ 

x cxp(-cr y ' oy oy ) x exp(--oy ay oy ) x exp(-— oy a y > oy ) x exp(-oy oy oy oy ), 
where the decomposition of each of the terms in Eq. 30 in elementary gates can be done using the methods de- 



scribed in previous works. 2,3 In particular, in Fig. 3 we show the decomposition of the term exp V*' 1 ^' 2 ^' 2 ^ 
where the qubits were relabeled as (n, i) = n + 2i — 1 (e.g., (0, 1) — ► 1). 

On the other hand, it is important to mention that the number of operations involved in the decomposition 
is not related to the distance between the sites i and j, as in the fermionic case. 



5. MEASUREMENT: CORRELATION FUNCTIONS AND ENERGY SPECTRA 

In previous work 2, 3 we introduced an efficient algorithm for the measurement of correlation functions in quantum 
systems. The idea is to make an indirect measurement, that is, we prepare an ancilla qubit (extra qubit) in a 
given initial state, then interact with the system whose properties one wants to measure, and finally we measure 
some observable of the ancilla to obtain information about the system. Particularly, we could be interested in 
the measurement of dynamical correlation functions of the form 



G(t) = (ip\T^A\TBjip) 



(31) 



where A t and Bj are unitary operators (any operator can be decomposed in a unitary operator basis as A = 
&iAi, B — YJ fijBj), T = e~ lHt is the time evolution operator of a time-independent Hamiltonian H, and 

i 3 

is the state of the system whose correlations one wants to determine. If we were interested in the evaluation of 
spatial correlation functions, we would replace the evolution operator T by the space translation operator. In 
Fig. 4 we show the quantum algorithm (quantum network) for the evaluation of G(t). As explained before, 2,3 
the initial state (ancilla plus system) has to be prepared in the quantum state |+) a ® (where a denotes 
the ancilla qubit and |+) = ^^2). Additionally, we have to perform an evolution (unitary operation) in the 
following three steps: i) a controlled evolution in the state |1) of the ancilla C-B = |0)(0| <g> I + |1)(1| <8> Bj, ii) 
a time evolution T, and iii) a controlled evolution in the state |0) of the ancilla C-A = |0)(0| ® A{ + ® 



Finally we measure the observable (2er+) = (tr^ + : 



G(t). 
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Figure 4. Quantum network for the evaluation of G(t) = (ip\T^ A^TBjip) . 



On the other hand, sometimes we are interested in obtaining the spectrum (eigenvalues) of a given observable 
Q (i.e., an hermitian operator). A quantum algorithm (network) for this purpose was also given in previous work. 3 



Again, the basic idea is to perform an indirect measurement using an extra qubit (see Fig. 5). Basically, we 
prepare the initial state (ancilla plus system) | +) a £g) | cj>) , then apply the evolution e 4< 3 <T ° i , and finally measure the 
observable (2a a + (t)) = (0|e -4< 3'<?!>). Since the initial state of the system can be written as a linear combination 

of eigenstates of Q, \<f>) — In |Vv»), where j n are complex coefficients and \ip n ) are eigenstates of Q with 

n=0 

eigenvalue A„, the classical Fourier transform applied to the function of time (2<7+(f)) gives us A„ 

L 

F{\) -^27r| 7 „| 2 (5(A-A„). 

n=0 



(32) 



Without loss of generality, we can choose Q — H, with H some particular Hamiltonian. 

It is important to note that in order to obtain the different eigenvalues of Q, the overlap between the initial 
state and the eigenstates of Q must be different from zero. One can use different mean-field solutions of Q as 
initial states \4>) depending on the part of the spectrum one wants to determine with higher accuracy. 
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Figure 5. Quantum network for the evaluation of the spectrum of an observable Q. 



6. ALGORITHM EFFICIENCY AND ERRORS 

An algorithm is considered efficient if the number of operations involved scales polynomially with the system 
size, and if the effort required to make the error e in the measurement of a relevant property smaller, scales 
polynomially with 1/e. 

While the evolution step involves a number of unitary operations that scales polynomially with the system 
size (such is the case for the Trotter approximation) whenever the Hamiltonian H is physical (e.g., is a sum of 
a number of terms that scales polynomially with the system size), the preparation of the initial state could be 
inefficient. Such inefficiency would arise, for example, if the state \ip) defined in Eq. 8 or Eq. 27 is a linear 
combination of an exponential number of states (L <~ x N , with N the number of sites in the system and x a 
positive number). However, if we assume that \ip) is a finite combination of states (L scales polynomially with 
N), its preparation can be done efficiently. (Any (Perelomov-Gilmore) generalized coherent state can be prepared 
in a number of steps that scales polynomially with the number of generators of the respective algebra.) On the 
other hand, the measurement process described in section 5 is always an efficient step, since it only involves the 
measurement of the spin of one qubit, despite the number of qubits or sites N of the quantum system. 

Errors e come from gate imperfections, the use of the Trotter approximation in the evolution operator, and 
the statistics in measuring the spin of the ancilla qubit (sections 3.2, 4.2, and 5). A precise description and study 
of the error sources can be found in previous work. 2 The result is that the algorithms described here, for the 
simulation of physical systems and processes, arc efficient if the preparation of the initial state is efficient, too. 



7. CONCLUSIONS 



We studied the implementation of quantum algorithms for the simulation of an arbitrary quantum physical 
system on a QC made of qubits, making a distinction between systems that are governed by Pauli's exclusion 
principle (fermions, hard-core bosons, anyons, spins, etc.), and systems that are not (e.g, canonical bosons). For 
the first class of quantum systems, we showed that a mapping between the corresponding algebra of operators and 
the spin-1/2 algebra exists, since both have a finite-dimensional representation. On the other hand, the operator 
representation of quantum systems that are not governed by an exclusion principle is infinite-dimensional, and an 
isomorphic mapping to the spin-1/2 algebra is not possible. However, one can work with a finite set of quantum 
states, setting a constraint, such as fixing the number of bosons in the system. Then, the representation of 
bosonic operators becomes finite-dimensional, and we showed that we can write down bosonic operators in the 
spin-1/2 language (Eq. 24), mapping bosonic states to spin-1/2 states (Eq. 23) . 

We also showed how to perform quantum simulations in a QC made of qubits (quantum networks), giving 
algorithms for the preparation of the initial state, the evolution, and the measurement of a relevant physical 
property, where in the most general case the unitary operations have to be approximated (sections 3.2,4.2). 

The mappings explained are efficient in the sense that we can perform them in a number of operations 
that scales polynomially with the system size. This implies that the evaluation of some correlation functions 
in quantums states that can be prepared efficiently is also efficient, showing an exponential speed-up of these 
algorithms with respect to their classical simulation. However, these mappings are insufficient to establish that 
quantum networks can simulate any physical problem efficiently. As we mentioned in the introduction, this is the 
case for the determination of the spectrum of the Hamiltonian in the two-dimensional Hubbard model, 3 where 
the signal-to-noise ratio decays exponentially with the system size. 

Finally, in Fig. 6 a table displays the advantages of simulating some known algorithms with a QC than with 
a CC, concluding that QCs behave as efficient devices for some quantum simulations. 



Algorithm 


Speed-Up • 


Correlation functions 


Exponential 


Factoring (Shor) 


Superpolynomial 


Search (Grover) 


Quadratic 



Figure 6. Quantum vs. classical simulations. Speed-up refers to the gain in speed of the quantum algorithms compared 
to the known classical ones. 
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